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Abstract. We prove that a potential q can be reconstructed from the 
Dirichlet-to-Neumann map for the Schrodinger operator — A g + q in a 
fixed admissible 3-dimensional Riemannian manifold (M,g). We also 
show that an admissible metric g in a fixed conformal class can be con- 
structed from the Dirichlet-to-Neumann map for A s . This is a construc- 
tive version of earlier uniqueness results by Dos Santos Ferreira et al. 
7 on admissible manifolds, and extends the reconstruction procedure 
of Nachman [21] in Euclidean space. The main points are the derivation 
of a boundary integral equation characterizing the boundary values of 
complex geometrical optics solutions, and the development of associated 
layer potentials adapted to a cylindrical geometry. 



1. Introduction 

This paper is concerned with the problem of reconstructing material pa- 
rameters of a medium from boundary measurements. A typical question of 
this type is Calderon's inverse conductivity problem [6], which consists in 
recovering the conductivity of a body from voltage to current measurements 
at the boundary. For bounded domains in Euclidean space in dimensions 
n > 3, it was proved in [30] that a smooth positive scalar conductivity 
a is uniquely determined by the Dirichlet-to-Neumann map (DN map) K a 
representing the boundary measurements. This uniqueness result was then 
extended to a reconstruction procedure in [21] and independently in [23], see 
also [13] . In two dimensions, uniqueness and reconstruction for this problem 
was proved even for bounded measurable conductivities in [2], [3]. 

In this paper we consider Calderon's inverse problem and related questions 
in anisotropic media, where the conductivity depends on direction. This 
corresponds to replacing the scalar conductivity a by a smooth symmetric 
positive definite matrix. The question then is to recover the matrix a from 
the DN map A CT , up to the natural obstruction given by diffeomorphisms 
which fix the boundary. If n = 2, it is proved in \\\ that any bounded 
measurable matrix conductivity a is determined by A CT up to diffeomorphism. 
For constructive results see [32] and the references therein. 

In three and higher dimensions the anisotropic Calderon problem is open 
even for smooth matrix conductivities. We refer to [7] for a more thorough 
discussion and references to known results. It was observed in [19] that the 
anisotropic Calderon problem is closely related to certain inverse problems 
for the Laplace-Beltrami operator on a Riemannian manifold, which we set 
out to define. 
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Statement of main results. Let (M, g) be a compact oriented Rie- 
mannian manifold with C°° boundary, and let A g be the Laplace-Beltrami 
operator. In local coordinates 

where g = (gjk) is the metric in local coordinates, (g^ k ) is the inverse matrix 
of (gjk), and \g\ = det(gjk). Consider the Dirichlet problem 

A g u =0 in M, 
u = f on dM. 

For any / G H 3 / 2 (dM) there is a unique solution u G H 2 (M), and the DN 
map is defined by 

A g : H 3 ' 2 (dM) H l l 2 (dM), f ^ d y u\ dM 

where the normal derivative is given by 

du 
dxj 

j,k=i J 

Here = Y^l 1=1 9ki u ' \ an d (v 1 , ■ ■ ■ , v n ) is the coordinate expression for the 
unit outer normal vector v on dM. 

Our first result states that the map A g constructively determines g within 
a known conformal class of admissible metrics (as defined below). 

Theorem 1.1. Let (M,g) be a given admissible 3-dimensional manifold. If 
c is a smooth positive function on M, then from the knowledge of A cg one 
can constructively determine c. 

The next question concerns an inverse problem for the Schrodinger equa- 
tion in (M,g). If q is a smooth function on M, we consider the Dirichlet 
problem 

(-A 9 + q)u = in M, 
u = f on dM. 
We make the standing assumption that 

is not a Dirichlet eigenvalue of — A g + q in M. 

This means that for any / G H 3 / 2 (dM) the equation has a unique solution 
u G H 2 (M), and the DN map can be defined by 

Ag, q : H 3 / 2 (8M) -> H l l 2 (dM), f ^ d v u\ aM . 

The second main result is as follows. 

Theorem 1.2. Let (M,g) be a given admissible 3-dimensional manifold. 
If q is a smooth function on M, then from the knowledge of A 9j(J one can 
constructively determine q. 

To complete the statement of the main results, let us give the definition 
of admissible manifolds. These arose in [7j as the first class of manifolds 
beyond real-analytic or Einstein ones for which one can prove uniqueness 
results for the anisotropic inverse problems described above. 
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Definition. A compact oriented Riemannian manifold (M, g) with smooth 
boundary is admissible if dim(M) > 3 and if (M,g) CC (T mt ,g) where 
T = R x Mq is a cylinder with metric g = c(e go) (here c is a smooth 
positive function and e is the Euclidean metric on R), and (Mo, go) is an 
(n — l)-dimensional simple manifold. 

Thus, up to a conformal factor, an admissible manifold is embedded in 
a cylinder (R x Mo,e © go) and therefore has a Euclidean direction. The 
transversal manifold (Mo, go) needs further to be simple: 

Definition. A compact manifold (Mo, go) with boundary is simple if for 
any p £ Mo the exponential map exp p is a diffeomorphism from its maximal 
domain in T p Mo onto Mo, and if dMo is strictly convex (meaning that the 
second fundamental form is positive definite). 

Examples of admissible manifolds include subdomains of the model spaces 
(Euclidean space, sphere minus a point, hyperbolic space), sufficently small 
subdomains of any conformally flat manifold, and domains in R n equipped 
with a metric of the form 



where c is a positive smooth function and where go is a simple metric in the 
x' variables. 

The unique determination results corresponding to Theorems 11.11 and 11.21 
were proved in [7J , where also earlier work is discussed and further references 
are given. There are several recent results that are concerned with the two- 
dimensional case. If M is a domain in R 2 , [5j proved the uniqueness result 
corresponding to Theorem 11.21 and briefly discussed a reconstruction pro- 
cedure. A proof with constructive character was given in |11| for arbitrary 
Riemann surfaces (M, g) with boundary (also for the magnetic Schrodinger 
operator), based on earlier nonconstructive proofs in [ID], [17]. Various 
constructive results for the two-dimensional case, also on Riemann surfaces, 
appear in [12], [H], [IS], [ID], and for the three-dimensional case an improved 
reconstruction result is given in [25J. 

Outline of proof. As mentioned, the uniqueness results corresponding 
to Theorems 11.11 and 11.21 were proved in [7J , but the proofs were not con- 
structive. The main point of the present paper is to give constructive proofs, 
following the well-known reconstruction procedure of Nachman |21] in the 
case where M is a bounded domain in R n and g is the Euclidean metric. 
We do not make any claims about the practicality of the reconstruction pro- 
cedure, but we do prove that all the steps in the corresponding uniqueness 
proofs in [7] can be carried out in a constructive way. 

The argument in |21] involves complex geometrical optics (CGO) solutions 
u = e~^' x (l + r) to the Schrodinger equation (— A J r q)u = in R n , and relies 
in a crucial way on a uniqueness notion for these solutions upon fixing decay 
at infinity. One has several equivalent ways of characterizing these solutions, 
and in particular it is possible to recover the boundary value u\qm as the 
unique solution to a boundary integral equation on dM involving A 9j(? and 
other known quantities. 
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CGO solutions on admissible Riemannian manifolds were constructed in 
[7J by Carleman estimates. The solutions were given in a compact manifold, 
and the construction did not involve a notion of uniqueness. The paper [18] 
introduced a direct Fourier analytic construction of CGO solutions, valid 
in the cylinder T and with a uniqueness notion obtained by fixing a decay 
condition in the Euclidean variable and Dirichlet boundary values on dT. 
We shall use the solutions constructed in [18 1 to prove Theorem 11.21 (which 
implies Theorem 1 1 . 1 1 aft er a simple reduction). 

We next sketch the proof of Theorem 11.21 Let (M,g) CC (T mt ,g) be an 
admissible manifold, and assume that g = e® go where (Mo, go) is simple. 
Here we suppose, for simplicity, that c = 1. We assume that A 9j(J and (M, g) 
are known (thus also A 9; o is known), and use the basic integral identity 



(Ag tq - Ag j0 )(u\ 9M )v dS = / quvdV (1.1) 
dM Jm 

which is valid for any solutions u,« £ H 2 (M) of (— A g + q)u = and 
— A g v = in M. We take u and v to be suitable CGO solutions such 
that u\qm may be obtained from A 9i(? as the unique solution of a boundary 
integral equation, and v\qm is explicitly given. Then the left hand side of 
(jl.ip is known. Taking the limit as r — > oo and varying certain parameters 
in the solutions, we recover in this way the quantity 



oo 

e -2Xr 



oo 

e 2iXxi q(x 1 ,r,9)dx 1 



dr (1.2) 



for any 9 £ S n 2 , where A is any nonzero real number and (r, 9) are polar 
normal coordinates in Mq with center on dMo- We have extended q into 
T \ M as a function in C c °°(T int ). 

Now, since (r, 9) are polar normal coordinates, the curves 7 : r 1— >■ (r, 9) 
are unit speed geodesies in (Mo, go) for any fixed 9. Denoting the quantity 
in brackets in (|1.2|) by f\(r,9), it follows that we have recovered 



e- 2Xr fx(l(r))dr 



„ -2 Ar 

'7 

for any maximal geodesic 7 going from dMo into Mo, and for any A / 0. 
This is the attenuated geodesic ray transform of f\ with constant attenuation 
— 2A, see [7], [29J. For any A such that this transform can be inverted, we 
recover f\ which is just the one-dimensional Fourier transform 

^ X1 {q(-,x' )}(-2A) forallx'GMo. 

It was proved in [7] and [8] that the attenuated ray transform is invertible 
for small |A|, thus giving information on ^ X i{q( • ; %')} for small frequencies. 
This determines the compactly supported function q( ■ ,x') uniquely by the 
Paley- Wiener theorem. To make this step more constructive one would like 
to invert the attenuated ray transform for all A, which would yield q( ■ , x') 
by taking the inverse Fourier transform. Up to now the argument has been 
valid for dim(M) > 3. However, if dim(M) = 3 then (Mo, go) is a 2D simple 
manifold, and the recent result [28] shows that the attenuated ray transform 
is invertible for any attenuation. We can use the inversion procedure in [28] 
to conclude the reconstruction of q from A 9j(? if (M, g) is 3-dimensional. 
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Boundary integral equation. The main new point in the proof is a 
Fredholm boundary integral equation characterizing the values on dM of 
suitable CGO solutions in T. This equation has the form 

(ld + 1 S T (A g>q -A gt0 ))f = u on dM. (1.3) 

Here uq is an explicit function on dM depending on various parameters, 
7 is the trace operator H 2 (M) — > H 3 / 2 (dM), and S T is a special single 
layer potential depending on a large parameter r and adapted to the CGO 
solutions and the geometry of the cylinder (T, g) = (M x Mo, e® go). In fact, 
we have 

S T f(x) = [ K T (x, y)f{y) dS(y), x G T int \ dM, 
JdM 

where the integral kernel K T (x,y) is explicitly determined by r and by the 
Dirichlet eigenvalues and eigenfunctions of the Laplace-Beltrami operator 
on (Mo, go)- We establish basic properties of the single layer operator in 
Section [2j 

We prove that for suitable choices of uq and for |r| large, the equation 
(|1.3p has a unique solution / G H 3 / 2 (dM), and one has / = u\qm where 
u is the corresponding CGO solution. Since the operator on the left hand 
side of (|1.3p is determined by the boundary measurements and since uq is 
explicit, we can indeed determine the boundary values of CGO solutions by 
solving this Fredholm integral equation. 

This approach is analogous to [21] which considers the Euclidean case, 
except that the uniqueness notion for CGO solutions is obtained from decay 
conditions and Dirichlet boundary values on the cylinder (T, g) instead of 
a decay condition at infinity in W 1 . Recently in [23J another constructive 
approach appeared. There the boundary integral equation is obtained via 
Carleman estimates in M, and no extension of M to a larger set is needed. 
It is presumable that a similar approach would work in our case. However, 
the single layer potential obtained from Carleman estimates is perhaps not 
so explicit as the operator S T introduced above. 

We remark here that it would be interesting to establish reconstruction 
results corresponding to Theorems 1 1 . 1 1 and 1 1 . 2 1 for the magnetic Schrodinger 
equation or for the time-harmonic Maxwell equations. Uniqueness results 
for these equations are proved in [7J, [IS]) and constructive results in the 
Euclidean case appear in [26], |27| . 

Structure of paper. Section [1] is the introduction. The basic proper- 
ties of the single layer operator S T and related Faddeev Green functions are 
considered in Section [2j In Section [3] we introduce several equivalent ways 
of characterizing CGO solutions, including the required boundary integral 
equation. The results in Sections [2] and [3] are in fact valid for any transver- 
sal manifold (Mo, go) with smooth boundary (not necessarily simple). The 
proofs of Theorems 11.11 and 11.21 are given in Section 01 

Acknowledgements. C.E.K. is partly supported by the NSF grant DMS- 
0968472. M.S. is supported in part by the Academy of Finland. G.U. is 
partly supported by NSF, a Chancellor Professorship at UC Berkeley and a 
Senior Clay Award. 



RECONSTRUCTIONS ON ADMISSIBLE MANIFOLDS 



6 



2. Boundary layer potentials 

Notation and function spaces. In this section we assume that (M, g) 
is a compact manifold with smooth boundary, having dimension n > 3, and 
that (M, g) CC (T int ,g) where T = R x M and g = e g , and (M , (fa) is 
any compact (n — l)-dimensional manifold with boundary (no restrictions 
on the metric go)- Points of T are written as x = {x\,x') where X\ is the 
Euclidean variable and x' is a point in the transversal manifold Mq. 

We write ( • , • ) for the inner product of tangent vectors, 1-forms, and 
other tensors, and | • | for the norm. The volume element in (T, g) is 

dV(x) = dV g {x) = dx! dV go (x') 

with dV go the volume element in (Mo,<7o)- We also write T = dM CC T, 
and denote by dS the volume element on T. 

Let L 2 (T) = L 2 (T, dV) be the standard I? space in T, and let H S (T) be 
the corresponding L 2 Sobolev spaces. Since Mo is compact, we define 

tffocCn ={/;/€ H s ([-R, R] x M ) for any R > 0}, 
H S C {T) = {/ G H S (T) ; f{x u x) = when |xi| > R for some R > 0}. 

Writing (t) = (1 + 1 2 ) 1 / 2 , we introduce for 6 / 0, s > the weighted spaces 
L 2 (T) = {feL 2 oc (T);(x 1 ) s feL 2 (T)}, 
H§(T) = {/ € #f oc (T) ; G ^ S (T)}, 

4o( T ) = {/e^(r);/|aT = 0}, 

^l 1 oc,o( T ) = U G ^locm;/l9T = 0}, 

and 

^-oo,o(r) = |J Hl (T). 
<5eR 

Also, Hq(T) = {/ G fl^T) ; f\ 9T = 0}. We define, in the L 2 (T) duality, 

On T = dM we consider the usual space L 2 (T) = L 2 (T,dS) and the corre- 
sponding Sobolev spaces H S (F). 

Let < Ai < A2 < • • • be the Dirichlet eigenvalues of — A go in (Mo, go), 
and let {<fii}fZi be an orthonormal basis of L 2 (Mo) consisting of Dirichlet 
eigenfunctions, 

-Agofa = W in M , faeH^Mo). 

We write Spec(— A 9o ) = {AJ^ 1 . If / G L 2 (T) we consider the partial 
Fourier coefficients 

f(xi,l) = (f(xi, - )Ai)l^{m )- 
The Parseval identity implies that 

/oo 00 
2l/>i,0l 2 ^i- 
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Standard single layer operator. Let Kq be the usual inverse of the 
Dirichlet Laplacian on T, defined as follows: since any u G Hq(T) satisfies 
the Poincare inequality 

p poo (• 

I \u\ 2 dV = I I \u(xi, x')\ 2 dV go dx\ 

JT J -co JM 

<[ c[ \d x ,u(x 1 ,x')\ 2 dV go dx 1 <C [ \du\ 2 dV, 

the bilinear form B(u, v) = f T (du, dv) dV is bounded and coercive on Hq(T). 
Consequently, for any F G if _1 (T) there is a unique weak solution u = KqF 
in Hq(T) of the equation — A g u = F in T. Thus Kq is a bounded linear 
operator 

K : H-\T) -> Hq(T), -A g K = Id. 
We also consider the trace operator which restricts functions to T, 

7: Hl(T) -> H 1 ' 2 ^), 1U = u\ r . 
The adjoint of 7 satisfies 

7* : ii" -1 / 2 (r) — >■ ij _1 (r), 

(l*h,(f) L 2 (T) = {h,j(p) L 2 {r) for any y> G #o( T )- 
Thus formally j*h = /idS". 

Definition. The standard single layer operator on T is the map 

So = K 7* :H- 1 l 2 {T)^Hl(T). 

The next result gives the basic jump and mapping properties of Sq. Here 
we write M_ = M , M + = T \ M, and 7^ = (u|jvf T )|r for the restriction 
of u to r from the interior or exterior. If u is a function on T \ T such that 
■u|m t are i? 1 in M T and satisfy — A g u = in M T , we define the normal 
derivatives from the interior or exterior weakly as elements of .f/" _1 / 2 (r) by 

({d u u) T ,h) L 2 ir) = ±{d{u\M T ),de h ) L 2 {MT) , h G H 1/2 (T), 

where G H 1 ^) is any function with e^|r = /1 and e/JdT = 0. The jumps 
on T are defined by 

[u] r = "i-u - 7+u, 
[d u u] r = {d v u)- - (<9„u)+. 

Lemma 2.1. If / G H~ l / 2 {T) then u = Sq/ G H%(T) satisfies 

-A g u = in M ± , 

Mr = 0, 
[^n] r = /• 

If / G H S (F) for s > -1/2 then u| M _ G H s+3 / 2 (M_) and u| M+ = u| M+ for 
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Proof. Harmonicity and the first jump property are direct consequences of 
the properties of Kq. The definitions also imply that for h G H l l 2 (Y), 

{[d v u\\ T ,h) L 2 {v) = {du,de h ) L2{T) = {7* f,e h ) L 2 {T) = (f,h) L2{r) . 

This shows the second jump property. If / G H k+1 / 2 (T) with k > 0, then u 
satisfies 

-A g ii = in M±, 

Mr = 0, 
[d„u] r = feH k + 1 / 2 (T). 

The transmission property [2Q|, Theorem 4.20] implies that u\m± are H k+2 
near T. Then the properties for u follow from standard interior and boundary 
regularity for — A g and from interpolation. □ 

Corollary 2.2. The trace single layer potential satisfies for s > —1/2 

7S0 :H S {T) ^H S+1 (T). 

Since Kq maps C£°(T mt ) to L 2 (T mt ), the Schwartz kernel theorem shows 
that there is a distributional integral kernel Ko(x,y). (We restrict to T int 
to avoid having to talk about distributions on manifolds with boundary.) 
Then formally 

S f(x) = J K (x,y)f(y)dS(y). 

We will need some basic properties of the integral kernel. These are easily 
obtained by comparing to the Green function on compact manifolds [3], |31j . 

Lemma 2.3. The kernel Kq is smooth in T mt x T mt away from the diagonal, 
and it satisfies (with d the Riemannian distance) 

\K Q (x,y)\ <Cud(x,y) 2 - n , x, y G U CC T int . 

Also, Ko(x,y) = Ko(y,x), and for any x £ T mt one has A g (Ko(x, •)) = 
in T int \ {x}. 

Proof. The condition Ko(x,y) = Kq(jj,x) follows since — A g is symmetric. 
Let U CC W CC T mt where U is open and (W,g) is a compact manifold 
with smooth boundary, and let G(x, y) be the Dirichlet Green function for 
the Laplacian in (W,g) [H Section 4.2]. We will prove that 

K (x,y) = G(x,y) + R(x,y), x,y£U 

where R G C°°(U x U). Since the function G(x, y) has the stated properties 
by [H Theorem 4.17], they also follow for Kq by using smoothness of R and 
simple arguments. 

Consider 9? G C£°(U) and let v = uq — ui where uq = Kqlp and 

u 1 (x)=l G(x,y)<p(y)dV(y). 
Jw 

Then — A g ui = <p in W with iii|avF = 0. It follows that v G H l (W), 
Iblli/^VF) — C||v , IIl 2 (i/)) an d A g v = in W. By elliptic regularity and 
Sobolev embedding, 



90)i 
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This may be rewritten as 

V k x R(x,yMy)dV(y) < C k \\<p\\ L2{u) . 

iw 

Consequently \\V^R(x, ■ ) 1 1 z^ 2 (t/) ^ Cfc uniformly over x G U . Now R(x, y) = 
R(y,x), so indeed R is smooth in U x U. □ 

r-dependent single layer potential. In |18j (see also Proposition 13.11 
below) it is shown that for any r G K with |r| > 1 and r 2 ^ Spec(— A 
and given 5 > 1/2, one has a bounded linear operator 

G T :L^T)^tf£ 5 (T)ntfI 5)0 (T) 

such that u = G T f is the unique solution in i?l oo0 (T) of the equation 
e TXl (-A g )e- TXl u = f in T, for any / G L 2 (T). 
If r is as above, we define another operator 

K T :L 2 c (T)^H? oc (T)nH{ ocfi (T), 

K T f = e- TX1 G T (e TXl f). 

From the properties of G T , it immediately follows that K T is an inverse for 
the Laplacian: one has 

-A g K T = Id on L 2 C {T). 

The map K T is a Faddeev type Green operator on T. It differs from the 
standard Green operator Kq by having exponential factors in its integral 
kernel. Also, K T - K maps L 2 C (T) into C°°(T) since -A g (K T - K )f = 
for all / G L 2 (T). This suggests the following result. We will do the proof 
with some care to ensure that the boundary dT does not pose a problem. 

Lemma 2.4. K T = Kq + R T where R T is an integral operator with kernel 
in C°°(T x T). 

Proof. We prove the result for r > 1, the case r < —1 being similar. The 
expression for G T in |18[ Proposition 4.1] shows that for any / G L 2 {T) with 
5 > 1/2 we have 

oo 

i=i 

with convergence in L 2 _ S (T), where for fj, G M \ {0} we define 

1 

If] — fjL 

This implies that for / G L 2 (T) one has 

oo 

K T f( Xl ,x') = - e~ TXl K + vxi t t-vxM t ' h ■ , 0)K*i)<MA 

i=i 

oo 



7>(t) = J?" 1 \ ——v{ri) \ it), v G y'(R). 



K f(x h x') = -^[^ v ^^_ v ^/(•,/)](x 1 )^(x , ) 



i=i 

with convergence in L 2 oc (T). The second identity follows by solving the 
equation — A g u = f in T, / G L 2 (T), u\qt = 0, by taking the Fourier 
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transform in x\ and expanding in terms of Dirichlet eigenfunctions in x' . 
Note that Kq is obtained from K T by setting r = formally. 
We compute for v £ L^(R) 

oo 

[e~ T{ - Xl ' yi) m T (xi - 2/1, v^) - moOi - Z/l, V^M^i) 

-oo 

where for a > 



m T (t,a) = y 



oo 

itr) 



oc 



e ltT 7 dry. 

(w/-(r + a))(M7-(r-a)) 



We claim that 



e rr m T (r,a) -m (r,a) = \ _i '_ or ^ °' (2.1) 

If this holds, then using that r / \/Al by the condition r 2 ^ Spec(— A 50 ) we 
obtain 

/OO 1 

R T (x,y)f(y)dV(y) 
where R T (x,y) is the integral kernel 

i? T (x,y) = - ^ -4= e-^^-^OOMl/O- 

/ /TV 2 ^ Ai 

Since the sum is finite, the kernel is smooth in T x T as required. 
It remains to prove (|2.ip . Note that for r£l and a > 0, 

I roo e i(j]+ir)r 

n ' 1 2tt (fa + *r)i - a) (fa + »r)» + a) 



T 



'TV 

where the last expression is a contour integral over the curve 7 T fa) = r] + ir 
for 7/ E (—00,00), and 

1 e irz 

F z) 



2tt (iz — a)(iz + a) 
It follows that 

e~ rr m T (r,a) —mo(r,a) = I / — / I F(z)dz. 

\-'7t ■'To/ 

We may interpret the right hand side as a limit of integrals over closed 
rectangular contours since F(±R + it) where < t < r decays as R — >■ 00. 
Now F{z) is analytic in C x {±ia} with simple poles at ±ia, so (|2.ip follows 
from the residue theorem. □ 

Since R t u\qt = for any u, the preceding result implies that 
K T :H~\T)^Hl^{T). 
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Definition. If |t| > 1 and r 2 ^ Spec(— A 9o ), we define the r-dependent 
single layer potential 

S T = K Tl * : H-V^T) ^ H} OCj0 (T). 

The basic properties of S T follow immediately from the previous results. 

Lemma 2.5. Let |r| > 1 and r 2 Spec(-A So ). If / G J H^ 1 / 2 (r) then 
u = S T f G Hy oc0 (T) satisfies 

-A g u = inM±, 

Mr = 0, 

[d v u\v = /• 

If / G H S (T) for s > -1/2 then u\ M . E if s+3/2 (Af_) and u| M+ = for 
some u G ifj^ 3//2 (T) n H^ oc0 (T). The trace single layer potential satisfies 

~fS T : # s (r) -> if s+1 (r), s > -i/2. 

One has 

S T f(x) = J K T (x,y)f(y)dS(y) 

where the kernel K T (x,y) is smooth off the diagonal in T int x T int and 
A fl (ir T (a;, •)) = in T int \ {x}. 

3. Boundary integral equation 



We now describe four equivalent problems for characterizing the CGO 
solutions u: a differential equation (DE), integral equation (IE), exterior 
problem (EP), and boundary integral equation (BE). In this section we as- 
sume that (M,g) CC (T int ,g) is a compact manifold with smooth boundary 
and T = R x Mq, g = e(Bgo, and (Mo, go) is any compact (n — l)-dimensional 
manifold with smooth boundary. We use the notations in Section [2j 

As a first step, we quote the basic existence and uniqueness result con- 
cerning G T from |X8(, Proposition 4.3 and subsequent remark]. 

Proposition 3.1. Let 5 > 1/2 and A / 0, and suppose that q G L°°(T) is 
compactly supported. There exists tq > 1 (if q = then to = 1) such that 
whenever 

|t| > to and t 2 ^ Spec(— A 9o ), 
then for any f = fi + f 2 where f x G L 2 {T), f 2 G L 2 _ S (T), and 

^x\Hi. ' -> x ) nas support in {|£i| > |A|} for a.e. x' G Mq, 
there is a unique solution w G Hz.^ q(T) of the equation 

e TXl (-A + q)e~ TXl w = f in T. 
Further, w G H X _ SQ (T) n H 2 5 (T), and one has tu = G r u where 
(Id + qG T )v = /, ||(Id + gGr)" 1 ^^ < 2. 
Finally, it; satisfies the estimates 



^6 



w 



H S _ S {T) < C\t\ S 1 ||/i||l2 (t) + ||/ 2 || L 2 ( T ) 



< s < 2, 



with C independent of r and /i, / 2 . 
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In the following result we extend q G L°°(M) by zero into T, and let to 
be as in Proposition 13.11 The result considers CGO solutions of the form 

u = u + e~ TX1 r 

where uq is any harmonic function in H? oc (T), and r G Hi (T). We will 
only fix the choice of the free solution uq in the next section. 



Proposition 3.2. Let q G L°°(M) be such that is not a Dirichlet eigen- 
value of — A g + q in M, and let \r\ > tq and r 2 ^ Spec(— A go ). Further, let 
u o £ H\ oc (T) be such that A g uo = in T. Consider the following problems: 

(DE) 



(IE) 



(EP) 



(BE) 



(- A s + q)u = in T 

e^(n-u )G J ffI oo>0 (T), 

u + K T (qu) = no in T 

i) A g u = in T \ M 

ii) u = u\t~^m for some u G Hf oc (T) 

iii) e TXl (n - Uo)|tnM = ^|r\M for some r G i^l^ oC 7 ") 

iv) (9^n) + = A 9j(? (7 + m) on T, 

(Id + 7S' T (A 9i9 - A 5 ,o))/ = i*o on T 

/e^ 3/2 (r). 



Each of these problems has a unique solution. Further, these problems are 
equivalent in the sense that u solves (DE) iff u solves (IE) , if u solves (DE) 
then u|t\M solves (EP), if u solves (EP) then there is a solution u of (DE) 
with G|t\M = u, if u solves (DE) then / = u\r solves (BE), and finally if / 
solves (BE) then there is a solution u of (DE) with u\r = f ■ 

Proof. The function u = uq + e~ TXl r solves (— A g + q)u = in T if and only 
if 

e TXl (-A g + q)e- TXl r = -e TXl qu in T. 

The right hand side is in L 2 (T), so by Proposition 13.11 there is a unique 
solution r G Hi^ (T). This proves that (DE) has a unique solution. It 
remains to prove that all four problems are equivalent in the sense described 
above. 

(DE) =^ (IE): Assume u solves (DE). Then u = uq + e~ TX1 r where 
r G fll^oCT), and 

e TX1 (-A g + q)e- TXl r = -e TXl qu in T. 

By Proposition 13. II we have r = G T v G H^ oc (T) where v satisfies 

v + qr = —e TXl qiiQ. 

Since q is compactly supported in T also v = —q(r + e TXl uo) is compactly 
supported. Thus we may apply G T to both sides of the last identity to 
obtain 

r + G T (qr) = -G T (e TX iqu ). 
Multiplying by e~ TX1 and adding uq to both sides gives (IE). 
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(IE) ==> (DE): Assume u solves (IE). Then the function r = e TX1 (u — uo) 
satisfies 

r = -G T (e TXl qu). (3.1) 

This shows that r G (T), and (DE) follows by applying — A g to both 
sides of (IE). 

(DE) ==>• (EP): Let u solve (DE), and define u = u\t\M- Clearly 
properties i), ii) and iii) of (EP) are valid. We need to show iv). Since u 
solves the equation (— A g + q)u = in M, we have 

(d v u) + = d u u\ T = A g<q (u\ r ) = Ag^+u). 

(EP) (DE): Suppose u solves (EP). Define v G H 2 (M) as the unique 
solution of the equation (— A g + q)v = in M with v\r = 7+w|r, and define 

v(x), x € M, 
u(x), x € T \ M. 

Then 7_n|r = 7+^|r and 

(d v u)-\r = A 3 , g (7 + u|r) = (d u u) + \ r 

by (EP) iv). It follows that u G H? oc (T) and (-A g + q)u = in T. Further, 
e r*i (u - u ) £ Fi^^CT) by (EP) iii). 

(DE) (BE): Let u solve (DE), and let / = u\ T . We fix a point 

x G T mt \ M and let v(y) = K T (x,y) where y G M. This is a smooth 
function in M by Lemma 12.51 

Now Green's theorem implies 

{ud v v — vd v u) dS = (uAgV — vA g u) dV. 
Jm 

By (DE) we have A g u = qu and d u u\r = A fl)9 /. Using the properties in 
Lemma 12.51 we obtain 



u(x) 



J ud u v dS — S T A gtq f(x) = —K T (qu)(s 



which is valid for x G T int \ M. The function v is harmonic in M, hence 
d v v\r = A g fi(v\r)- The symmetry of A g> o implies 

J ud v v dS = j uA 9; o(^|r) dS = j A g ^(u\r)v dS = S T A g) of(x). 

We obtain 

S T {A m -A gfl )f = K T {qu) in T int \M. (3.2) 

Adding u to both sides, using the fact that u solves (IE), and taking traces 
on r gives (BE). 

(BE) =>- (EP): Let / solve (BE). We define a function u G fl£ c (T) by 

u = u - S T (Ag jq - A Sj0 )/. 

This function is harmonic in T \ T by Lemma 12.51 and u\r = f by using 
(BE). The jump relation for S T implies that on T 

{d v u)- - (d u u) + = -(A g>q - A g ,o)f. 

But (d u u)- = A g fif, so we have (d u u) + = A 5j q(7 + tt). Therefore u\t^m 
satisfies (EP) i) and iv). Also (EP) ii) is valid by mapping properties of S T . 
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To prove (EP) iii) it is sufficient to show that for any h G ff 1//2 (r), 

e TXl S T h\T^M = w\t^m forsomeweffl^^T). (3.3) 

Formally one has e TX1 S T h = G T e TX1 ^*h where G T maps L 2 (T) to (T). 
However, we have not proved that G T has good mapping properties on neg- 
ative order Sobolev spaces. Thus, the proof proceeds differently and in- 
volves an extension w of e TXl S t }i\t-^m into T such that w = G T ip for some 
ip G L 2 {T). This will imply (|3.3|) by the mapping properties of G r . 
Define 

" e TXl (P (-fS T h) + F), x G M, 
e TX1 S T h, x£T\M 



w(x) 



where Pq : H 3 ^ 2 (T) — > H 2 (M) is the Poisson operator mapping ho to the 
function vq with —A g vo = in M and i>o|r = ^O) an d F G H 2 (M) is a 
function chosen so that e~ TXl w G H 2 oc (T). Clearly we need that F\r = 0, 
and since 

[d u (e- TXl w)} r = A g ,ohS T h) + d v F\ r - (d u S T h) + 

we also require that d u F\r = (d u S T h) + — A 9i o(7<SV/i) G H l / 2 (T). We can 
take F to be any function in H 2 (M) with this Cauchy data, and then e~ TXl w 
and also w is in H 2 oc (T). 
We now observe that 

K 3) 10, x G T \ M. 

Since w; G H 2 oc (T) this implies that e™ 1 (-A 9 )e" ra:i u; = ^ where ip G LcCO- 
Consequently ui = G T ip G -ffioo o 0O an d we have proved (|3.3p . □ 

Finally, let us verify that the boundary integral equation (BE) in Propo- 
sition is indeed Fredholm. 



Proposition 3.3. The operator 

7-SrCA,,, - a 9 ,o) : # 3/2 (r) -> # 3 / 2 (r) 

is compact. 

Proof. Let / G H 3 / 2 (T), and let u = P q f where P q : H 3 / 2 (T) -> F 2 (M) 
is the Poisson operator mapping ho to -uo where (— A g + g)-uo = in M 
and voir = ^o- The exact same argument leading to (|3.2p in the proof of 
Proposition 13.21 shows that 

S T (k m -k gfi )f = K T (qEJu) in T int \M 

where E : L 2 {M) -> L 2 (T) is extension by zero and J : H 2 (M) -)■ L 2 (M) is 
the natural inclusion. Taking traces on T, we obtain the factorization 

75 T (A 9i9 - A Sj0 ) = -fK T qEJP q . 

The result follows since on the right hand side J is compact and all other 
operators are bounded. □ 
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4. Proofs of the main results 

In Sections[2]and[3]we considered layer potentials and equivalent problems 
characterizing CGO solutions in the case where (M,g) CC (T mt ,g) where 
T = RxMo, g = e®go, and (Mo, go) can be any compact (n— l)-dimensional 
manifold with boundary. Now we specialize to the case where (Mo, go) is 
simple and prove Theorems 11.11 and 11.21 

The first step is to fix the harmonic functions uq used in Proposition 13.21 
We first choose a simple manifold (Mo, go) such that (Mo, go) CC (Mo, go). 
Below, we will write (r, 0) for the polar normal coordinates in (Mo, go) with 
center at a given point p 6 Mo \ Mo (these exist globally because the 
manifold is simple), and we write, following Section 5], 

a = a{x u r,0) = e - iTr \g\- 1/4 e iX{xi+ir) b(e) 

where A is a fixed nonzero real number and b £ C°°{S n ~ 2 ) is a fixed function. 
Note that a G C°°(T) since the coordinates (r, 9) are smooth in Mo. 

Proposition 4.1. Given any r with |r| > 1 and r 2 ^ Spec(— A go ), and for 
any point p G Mo x Mo, for any real number A / 0, and for any smooth 
function b = b(9), there is a function uq with 

A g u = in T, u € Hf oc (T), 

of the form 

u = e~ TX1 a + e~ TXl r 

where tq = G T f for some explicit function / and ||?"o||l 2 (a/) = 0(|r| _1 ) as 
|t| — > oo. 

Proof. If uq is of the required form, then A g uo = is equivalent with 

e™ 1 (-A 9 )e-™ 1 r = / (4.1) 
where / = e TXl A g (e~ TXl a). Writing $ = x\ + ir, we compute 

/ = - e - lTr [e T ' s '(-A g )e- T,1 '}(\g\- 1 / i e iXixi+ir) b(e)) 
= _ e -iTr[_ T 2^ $) ^ + r ( 2 (d$, d • ) + A 9 $) - A 9 ](|gr 1 / 4 e iA(xi+ir) 6(^)). 

Here we have extended ( • , • ) as a complex bilinear form to complex val- 
ued 1-forms. As in [TJ Section 5], we see that (d&,d$>) = and also that 
(2(d$,d-) + A 9 $)(|g|" 1 / 4 e iA ( :ri+ir )6(0)) = (this was the reason for the 
choice of a). Consequently 

/ = e- iTr A g (e irr a) = e~ irT (d\ + A go )(e irr a) 

= e^[e-^ r (A go - A 2 )(|g|- 1 / 4 e- Ar 6(0))]. 

Then for any 5 > 1/2 one has / G L 2 5 (T), the norm ||/||j^2 ^ is indepen- 
dent of r, and the Fourier transform ^ X if( ' , is supported in {|£i| > |A|}. 
By Proposition 13.11 we have a solution ro = G T f of (|4.ip . which gives the 
required solution no- □ 

We can now prove the main theorems. 
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Proof of Theorem We first consider the case, as in the beginning of this 
section, where (M, g) is an admissible manifold with conformal factor c = 1. 
Suppose that the manifold (M,g), and consequently also (Mo, go), and the 
map Ag jq are known. We wish to determine q from this knowledge. 
First note the basic integral identity (see [7J Lemma 6.1]) 

/ ((A gtg -A gfi )f 1 )f 2 dS= [ qu lU2 dV (4.2) 

JdM JM 

which is valid for any uj £ H 2 (M) with (—A + q)u\ = in M, A g u 2 = in 
M, and u^\qm = fj- We consider CGO solutions in T of the form 

ui = ii ,i + e" ra:i ri, 
U 2 = uo,2 

where uoj are harmonic functions provided by Proposition 14.11 having the 
form 

U01 = e - T(xi+ir) \g\- 1/4 e lX ^ +ir U(6) + e-^G^!, 
n0j2 = e T(*i+ir)| ff |-l/4 e a(z 1+ ir) +e ™i G _ r ^ 2 . 

Here r > To and r 2 ^ Spec(— A go ), (r, 6) are polar normal coordinates in 
(Mo, go) with center at p £ Md \ Md, A ^ 0, b is a smooth function in S n ~ 2 , 
and V'j are explicit functions with [ | f^rfcT- V'j 1 1 ^ (Af ) = 0(|r| _1 ). 

The point is that uoj are explicit functions which can be constructed from 
the knowledge of (M,g), and also f 2 = uo, 2 \r is known. By Proposition 13.21 
there is a unique CGO solution u\ of the above form, and the boundary 
value /i = u\\y is the unique solution in H 3 / 2 (T) of the boundary integral 
equation 

(Id + jS T (Ag t g - A Sj0 ))/i = «o,i on T. 

Since the operator on the left and the function on the right are known from 
our data, we can construct f\ as the unique solution of this Fredholm integral 
equation. Then the left hand side of fj4.2|) is known, and consequently we 
can determine from our data the integrals 

qu\u 2 dV (4-3) 

for any u\ and u 2 as above. 

Since u\ solves (— A g + q)u\ = with r\ G H^^^T), Proposition 13.11 
shows that 

r 1 = -G T (Id + (/ G T )- 1 ( e ™ 1 ^o,i) 
and using the form of iio,i and norm estimates for G T gives 

ih!l2(m) = ^(kr 1 )- 

Thus, taking the limit as r — > oo in (|4.3p . we have recovered from our 
boundary data the quantities 

[ q\g\- 1 / 2 e 2iX( - Xl+ir ^b(6)dV. 

JM 



\ 
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At this point it is convenient to extend q into T as a function in C^°(T mt ). 
This may be done by recovering the Taylor series of q on dM via bound- 
ary determination Section 8] (this procedure is constructive), and by 
extending q to a function in C£°(T mt ) so that £?|r\Af is known. Using that 
dV = \g\ 1 / 2 dx\dr dO, the last integral becomes 



e -2Ar 



oo 

e 2iXxi q(x 1 ,r,9)dx 1 



b(6)dr dO. 



10 JS"- 2 

Denoting the quantity in brackets by f\(r,0) and by varying the smooth 
function b, we determine the integrals 

roc 

2Ar r r n\ 7 r — 11 a r- an-2 



poo 

/ e~ 2Xr f x (r,9)dr for all 9 £ S 
Jo 



These integrals are known for any nonzero real number A and for any 
point p £ Mq \ Mq which is the center of the polar normal coordinates (r, 6) 
in Mq. Noting that r i-> (r, 6) is the unit speed geodesic in (Mo, go) starting 
at p in direction 9, and letting p approach 8Mq, we can recover the integrals 

cT 

e- 2Xr hin{r))dr (4.4) 

for any geodesic 7 : [0, T] — > Mq where 7(0), 7 (T) £ 8Mq and ~f(t) for 
< t < T lies in M nt . This is the attenuated geodesic ray transform of f\ 
in (Mo, go), with constant attenuation — 2A. Now, assuming dim(M) = 3 so 
(Mo, go) is 2-dimensional, we invoke the invertibility result for the attenuated 
ray transform [28] which allows to recover the function fx in Mo from the 
integrals (|4.4p for any A. Thus, we have determined the integrals 

00 

e 2iXxi q(x 1 ,x')dx 1 

-00 

for any A / and for any x' £ Mo. This determines q in M by inverting 
the one-dimensional Fourier transform. 

We have proved the theorem in the case where (M,g) is an admissible 
manifold and with conformal factor c = 1. For general conformal factors, 
suppose that (M, g) is admissible and g = eg where g = e (B go- Define also 
q = c(q — q c ) where q c = c~Aq(c ~). The identity 

n+2 n— 2 

c 4 (-A C g + q)(c 4 u ) = (-Ag + q)u 



implies that, since v„ = c 1//2 



< v g, 

z2l „n . n — 2 
T 

Thus, from the knowledge of A g ^ q and (M,g) we can determine Agq. The 
proof above then shows that one can reconstruct q, from which q is easily 
determined. □ 

Proof of Theorem li.il Let (M, g) be admissible and known and suppose 
that A C g = A C g t o is known. By boundary determination [7] we can determine 
c\qm and d u -c\gM- The identity (|4.5p shows that 

n n-2 Tl — 2 -1 

A M f = c*A C g t o(c 4 f) + —^ c 1 (d u ~c)f 



A M f = cUg^c-^f) + ^-j-c-i^c)/. (4.5) 
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with q = — c 4 A C g(c ± ). This shows that A cg determines Agg, and 
Theorem 11.21 implies that we can recover q. 
We write w = log c * and compute 

n 
n 

i,fc=i 

n 

n+2 n — 2 \ — ^ n — 2 - r, n — 2 n — 2 

= c— A cg (c- — )+ 2^ c— g^djic- — )d k (c —). 
j,k=i 

This implies that 

— A g w + (dw, dw) g = q in M, 

w|aM = log c 4 | 9M . 

This nonlinear Dirichlet problem has a unique solution by the maximum 

principle [U] and we have already recovered the right hand side q and the 

— n ~ 2 

boundary value log c 4 \g M , so we may construct w in M by solving the 
problem. This determines c in M. □ 
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